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SYNOPSIS 


Recent  investigations  on  groundshock  phenomena  indicate  that  the 
dynamic  response  of  certain  non* linear  (locking)  media,  may  approximate 
the  behavior  of  cohcoive  granular  soils  and  porous  rock  under  high  pres¬ 
sures.  There  is  a  considerable  body  of  fairly  recent  American  and  USSR 
scientific  literature  on  this  topic.  It  seems  that  the  Russian  investi¬ 
gations  parallel  our  efforts  in  this  field.  The  theory  of  wave  propaga¬ 
tion  in  such  materials  should  have  application  to  the  design  of  underground 
shelters  at  very  high  preooure  levels  and  may  shed  some  additional  light, 
by  these  relatively  simple  means  on  problems  associated  vlth  phenomena 
near  ground  zero,  in  case  of  surface  bursts  or  underground  explosions. 

In  Sections  II  and  III  results  of  current  investigations  are  presented. 


I.  RECENT  AND  CURRENT  INVESTIGATIONS. 


As  *  result  of  our  Investigations  into  groundshoek  phenomena  in  the 
recent  years,  ve  have  conjectured  that  the  behavior  of  porous  or  granular 
solids  (such  as  dry  sand  or  porous  rock)  could  he  characterized  by  a  med¬ 
ium  which  under  uniaxial  compression  exhibits  a  stress-strain  diagram  which 
tends  to  become  parallel  with  the  stress  axis  beyond  a  certain  critical 
strain,  (Figure  l).  Such  a  stress-strain  diagram  may  be  replaced  by  vari¬ 
ous  bi- linear  approximations,  and  such  materials  are  currently  referred 
to  as  locking  media,  or  materials  of  limited  compressibility,  (Figure  2a, 
b,  c,  d). 

The  static  behavior  of  such  materials  has  been  explored  to  a  consi¬ 
derable  extent  since  about  1957  (References  1,  2,  3,  4], 

The  dynamic  problem  for  the  case  of  the  upwards  curving  Btress- strain 
diagram  has  been  considered  previously,  [References  5  and  6],  but  the  assump¬ 
tion  of  a  stress-strain  curve  having  a  continuous  derivative,  (Figure  l) 
makes  it  difficult  to  obtain  closed  form  solutions.  Some  more  recent  ef¬ 
forts  in  this  direction  provide  additional  results  for  specific  initial 
conditions,  [Reference  7l» 

The  linear  approximations  shown  on  Figure  2a,  b,  c  have  resulted  in 
relatively  simple  solutions  and  some  of  these  results  have  been  applied 
with  some  success  to  the  prediction  of  groundshoek  phenomena,  [References 
8  and  9)  and  have  also  been  published  in  the  open  literature,  [References 
10  and  11] .  A  medium  which  shows  a  residual  elasticity,  (Figure  2d),  is 
treated  in  Section  II  of  this  report. 
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If  tbs  tvo  or  three  dimensional  ease  1>  considered,  one  may  replace 
the  stress-strain  diagram  of  the  uniaxial  oase  vith  a  stress-dilatation 
diagram  to  obtain  similar  results.  After  the  critical  ralue  of  the  dlla- 
tatlca  Is  reached,  a  post  compaction  yield  condition  must  be  specified. 

One  may,  for  example,  assume  that  the  medium  behaves  libs  an  invlscid 
fluid  or  that  it  follow  Coulombs  law  of  failure.  Buoh  a  medium  for  the 
ease  of  spherical  symmetry  1*  treated  in  lection  XXX  of  this  report. 

During  our  researches  into  this  sub Jest  ve  also  ooted  that  a  '  cnslder- 
able  interest  in  these  topios  exists  in  the  BBflR,  as  shown  by  Reference 
12,  13,  lk,  13,  16  and  17.  It  appears,  that  the  Russian  Investigations 
parallel  our  own  efforts,  and  the  extensive  numerical  work  of  Reference  17, 
indicates  that  the  application  of  these  researches  Is  probably  also  seri¬ 
ously  considered.  Significantly,  most  of  thoae  pepere  are  directed  towards 
problems  associated  vith  explosions  and  apparently  the  USSR  researchers 
suffer  from  lack  of  experimental  data  just  as  much  as  our  own  investiga¬ 
tions  do. 

Ae  the  result  of  our  own  work  and  the  review  of  the  USSR  literature. 

It  la  believed  that  these  Investigations  have  application  to  the  following 
topics  of  current  interest! 

(a)  Groundshook  phenemona  In  dry  granular  and  cohesive  soils,  and 
porous  rock. 

(b)  Feasibility  studies  for  the  design  of  underground  shelters  near 
ground  tero. 

(o)  Groundsbock  phenomena  in  the  neighborhood  of  ground  sero  of 
surface  bursts. 

(d)  Phenomena  associated  with  nuclear  explosions  contained  in 


underground  eavltlea 
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Nmnerical  evaluations  lc  a  fev  specific  cases  leads  us  to  believe 
that  this  approach  may  Indeed  contribute  to  the  solution  of  the  above 
problem,  [References  Id  and  19 ].  Zt  appears  that  the  Intrinsic  dissipative 
behavior  of  the  medium  is  the  crucial  characteristic  which  makes  it  approp¬ 
riate  for  such  investigations.  Further  the  tret  leal  and  experimental  work 
in  this  field  may  be  quite  fruitful. 

In  the  following  Sections  IZ  and  ZZZ  of  this  report  some  additional 
contributions  to  the  theory  of  locking  media  are  presented. 

In  Section  II  a  medium  with  a  stress-strain  diagram  as  shown  on  Fig¬ 
ure  2d  is  considered.  The  purpose  of  this  Investigation  is  to  clarify 
the  effects  of  "residual  elasticity”,  that  is,  elastic  action  which  persists 
after  compaction  or  locking  has  occurred.  This  modification  of  the  stress- 
strain  relationship  is  of  importance  since  it  brings  the  idealized  theory 
closer  to  the  anticipated  behavior  of  real  materials.  In  our  first  and 
second  report  to  the  MITRE  Corporation  [References  18  and  19 ]  the  effects 
of  this  behavior  were  numerically  evaluated,  and  it  was  shown  that  small 
residual  elasticity  after  compaction  does  not  significantly  influence 
the  attenuation  of  peak  stresses  in  the  medium  at  depth  of  500  feet  or 
more,  as  compared  to  a  locking  medium  without  such  residual  elasticity. 

This  medium  also  permits  us  to  draw  certain  conclusions  regarding  non¬ 
linear  effects  due  to  the  Influence  of  the  high  Intensity  airslap  in  the 
vicinity  of  ground  zero. 

In  Section  III  spherical  waves  in  an  ideal  locking  medium  are  in¬ 
vestigated.  These  results  prove  to  be  useful  in  clarifying  the  close-in 
effects  of  the  exploding  casing  at  the  point  of  impact.  Numerical  evalua¬ 
tions  of  this  work  will  appear  in  our  subsequent  reports  on  this  subject. 
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SECTION  II 


WAVE  PROPAGATION  IN  A  PLASTIC - ELASTIC  MEDIUM 


By 


Richard  Skalak 


I.  STATEMENT  OF  PROBLEM. 

The  response  of  a  seml-inflnlte  body  of  plastic-elastic  medium.  Fig¬ 
ure  1,  is  to  be  found  due  to  a  blast  vave  pressure  applied  to  the  surface 
of  the  medium.  The  stress-strain  curve  of  the  material  is  shown  In  Fig¬ 
ure  2.  The  pressure  loading  is  assumed  to  be  applied  to  the  entire  Bur- 
face  of  the  medium  so  that  the  problem  is  one-dimensional.  The  type  of 
time  history  of  the  surface  pressure  considered  is  shown  in  Figure  3.  The 
strains  are  assumed  to  be  small  throughout. 


II.  GOVERNINO  EQUATIONS. 

(a)  The  Shock  Wave. 

The  applied  loading  generates  a  plane  shock  wave  which  compress^ » 
the  medium  from  its  initial  state  to  a  state  corresponding  to  some  point 
on  the  inclined  portion  of  the  stress-strain  curve.  The  velocity  of 
propagation  of  this  shock  vave  at  any  time  is 
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where  E  is  the  secant  modulus  of  elasticity.  Figure  2 
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where  0  ,  c8  are  the  normal  vertical  stress  and  strain  immediately  behind 

the  shock  front,  the  density  p  is  assumed  constant  because  the  strains  are 
assumed  to  be  small. 

The  particle  velocity  behind  the  shock  wave  is 
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Za  terms  of  the  modulus  I  of  tho  inclined  portion  of  the  •trees- 
•train  curve,  Figure  2,  the  value  of  «B  can  be  vritten  ae 
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vhere  «e  is  the  plastic  strain  shown  in  figure  2,  and  og  -  ll/p. 
From  liquations  (l)  through  (4)  it  follows  thatt 


a  m  Of  s 
ss  e  s 


i  -  <r>* 


(5) 


(6) 


(b)  Region  Behind  the  Shock  Wave. 

Behind  the  shod  wave  the  state  of  the  medium  Is  assumed  to  cor  re*, 
•pond  to  points  on  the  inclined  portion  of  the  stress-strain  curve.  Fig¬ 
ure  2.  If  u  is  the  downward  displacement  of  eny  particle  from  its  orig¬ 
inal  position,  the  equation  of  motion  1st 


dffs  3u 

51""  “ 


By  definition,  considering  coepresslve  strains  positive, 


du 

«"  "3s 

Since  se  is  a  constant,  (8)  yields 


8(s  •  «e) 
dt 
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Frcoi  Equation  (9)  and 
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It  follow* 
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The  u  in  Equations  (7)  and  (11)  1*  th*  parti  cl*  velocity.  Froo  thesa 
equation a, 


2 
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III.  80LUTI0M  FOR  A  PIVOT  SURFACE  PRESSURE  HISTORY. 

Consider  a  general  field  point  (s,  t)  in  Figure  k  and  two  points 
labeled  1  and  2  on  the  shock  curve  located  by  the  characteristics  passing 
through  the  field  point  as  shown  in  Figure  4.  By  virtue  of  Equations  (7), 
(11)  and  (12),  the  entire  field  behind  the  shock  wave  has  the  fora 

at  *  f(*  "  ct)  +  s(*  +  ct)  (13) 

u  -  ^  £f(t  -  ct)  -  g(*  +  ct  )]  (1*0 

where  f  and  g  are  unknown  functions. 

Writing  Equations  (13)  and  (14)  for  the  point  2  in  Figure  4t 


V  ■  f<‘2 

-  ct2)  +  g(*2  ♦  ctg) 

(15) 

"a  •  1  [fl 

>2  "  ct2>  ‘  «(*2  ♦  ctgjj 

(16) 

rroa  (1$)  and  (1 6) 


f(*a  "  ^  £  ®«2  +  2c"  *2  (17) 

and  tinea  alone  thn  characteristic  (a  -  et)  la  constant,  It  follows  that 
for  the  field  point  a,  t  indicated  in  figure  k 

f(a  -  ct)  •  j  e^g  ♦  Ug  (10) 

Similarly,  considering  point  1  it  is  found  that  for  the  field  point  a,  t 

^(b  +  ct)  ■  j  \  (^9) 

The  solution  for  any  field  point  can  now  be  expressed  in  terms  of 
values  a  and  u  directly  behind  the  shock  front  by  substituting  (18)  and 


(19)  into  (13)  and  (14) i 

°tm\  <ffal  +  ®«2>  “  s§  (41  ‘  V 
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(21) 

The  location  of  the  shock  front  a  (t)  and  the  values  of  stress  and 
velocity  o  and  Cl  at  the  shock  front  are  unknown  in  advance  for  a  given 

Z0  B 

surface  pressure  history.  However,  if  the  location  of  the  shock  vave  a  (t) 
is  determined,  then  a  and  u  are  given  by  Equations  (5)  and  (6)  and  the 
remainder  of  the  field  can  be  found  from  (20)  and  (21).  The  solution  is 
completed  by  establishing  s  (t)  for  a  given  surface  pressure  history  as 
follows: 

Consider  a  field  point  which  is  on  the  surface  *  ■  0  as  shown  in 
figure  (5)»  for  this  point  the  stress  o.  given  by  Equation  (20)  is  the 
applied  pressure  p  at  tias  t.  By  substituting  the  values  for  o  and  u 

ZB  8 
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from  Equations  (5)  and  (6)  Into  Equation  (20),  an  equation  relating  p,  z^  , 
and  Kg  la  obtained,  which  after  slnpllfi cation  reads t 


ET  P 
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(22) 


If  the  point  (0,  t)  in  Figure  (?)  la  assumed  to  move  along  the  t  axis, 
the  pressure  p(t)  is  the  given  surface  pressure  history.  The  heights 
shown  In  Figure  (5)  labeled  t^  and  may  then  be  also  considered  to  be 
functions  of  t.  Similarly  the  times  t^  and  tg  may  be  considered  as  func¬ 
tions  of  t«  The  time  t^  Is  related  to  t  by 


*i  -  t  -  r 
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Differentiating  with  respect  to  t, 
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dt 

dz 
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The 


1 

derivative  is  what  is  meant  by  t^  in  Equation  (22),  for  example. 


From  Equation  (24) 


dt 
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Defining  z|  to  be  the  derivative  ,  it  follows,  using  Equation  (25), 


dt. 
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dz^  dt^ 
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Similarly,  using 
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it  is  found  that 
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From  Equations  (26)  and  (28) 


c 


c  +  x. 


Substituting  Equations  (29)  into  (22) 


(29) 


g;  »  *  H  -  4  <*» 

vhera  p  is  p(t)  and  the  primes  denote  differentiation  vlth  respect  to  this 
same  t.  Equation  (30)  can  then  be  Integrated  with  respect  to  t  to  yield 
the  final  equation 


ir  «*)  ■  *2  -  *i  (3D 

c 

where  l(t)  is  the  impulse  of  the  applied  pressure  per  unit  area  defined  by 

t 

I(t)  m  J  p  dt  (32) 

o 

The  boundary  condition  tg  ■  ^  ■  0  at  t  ■  0  was  also  used  in  deriving  (31). 

Equation  (31)  defines  the  locus  x  (t)  of  the  shock  location  in  the 
x,  t  plane  and  completes  the  solution  of  the  problem. 


IV.  COMPUTATIONAL  PROCEDURE. 

The  Equation  (31)  may  be  used  to  solve  for  the  location  of  the  shock 
fropt  x2  at  a  time  tg  when  the  position  x^  is  known  for  the  earlier  time 
t^.  To  begin  the  computations  using  Equation  (31)  it  is  necessary  to 
establish  x^  at  seme  early  time  t^  as  a  starting  point.  A  first  and  a 
second  approximation  are  given  below  and  higher  order  approximations  can 
be  developed  by  similar  procedures. 
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a)  first  Approximation  for  Short  Times. 


A  pressure  time  history  such  as  shorn  In  Figure  3  could  be  replaced 
by  a  constant  pressure  pa  for  a  short  tils  tQ  so  that  the  impulse  during 
this  Interval  Is  conserved.  Then  for  the  time  interval  selected  Is  a 
constant  and  may  be  confuted  from  Equation  (5) t 


“20 


(33) 


The  value  of  s£  Is  also  a  constant  s^  and  is  given  by  Equation  (28) .  Then 
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A  starting  point  is  now  established  by  using 


*l(t)  "  *20 


(35) 


and 


t  *  to  +  "c  (36) 

b)  Second  Approximation  for  Short  Times. 

A  pressure  time  history  such  as  shown  in  Figure  3  can  be  represented 
by  a  Taylor  series: 


t2 

p(t)  -  P0  +  P^t  +  p£  T  . .  (37) 

The  values  of  ^  and  may  also  be  written  i-n  series  form: 

*1  “  *1*  +  k2t2  +  *3^  +  •  .  •  «  (38) 

*2  "  “l*  +  “2t2  +  “3^  * 
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The  Impulse  applied  at  the  surface,  Equation  (32)  Is: 

♦2  *3 

l(t)  -  P0*  +  Pq  T  *  p"  +  •  •  •  •  (40) 

Substituting  (38),  (39)  and  (40)  Into  Equation  (30)  and  equating  coeffici¬ 
ents  of  each  power  of  t  yields: 


■1  *  *1 


(U) 


c 


etc 


The  functions  z^(t),  Equation  (38),  muBt  sl.so  be  such  that  at  the  appropri¬ 
ate  time  It  Is  equal  to  the  value  of  z g  of  an  earlier  time.  This  re¬ 
quirement  may  be  written  as: 

2z.(t) 

zL(t)  -  *g(t  -  -f—  )  (43) 

where  all  brackets  mean  "a  function  of".  Substituting  (38)  and  (39)  In 
Equation  (43)  yields 

+  kgt2  +  k3t3  +  .  .  .  .  -  m1(t  -  |(kjt  +  kgt2  +  k3t3  +...)) 

+  Bg(t  -  |(kxt  +  kgt2  +  k3t3  +...))2 

+  m3(t  "  fftjt  +  4gt2  +  k3t3  +...))3 

+  etc...  (44) 


Equating  coefficients  of  each  power  of  t  in  Equation  (1*1*),  after  expand¬ 
ing  each  term  gives 


2k. 

V1  •  -r> 

(1*5) 

2k-  2k.  g 

“  “l  —  4  “2(1  '  "c_) 

(1*6) 

etc... 

Equations  (fcl)  and  (1*5)  may  be  solved  for  1^,  b^j 


Equations  (1*2)  and  (1*6)  nay  be  solved  for  kg,  mg  considering  k^  known: 


“2 


cp'o 

k3 

*1 

“i  *  ki 

cp; 

4 

*c 

®i  *  ki 

(W 

(50) 


Equation  (39)  con  now  be  used  with  the  coefficients  m^,  mg  given  by  (1*8) 
and  (50)  to  derive  the  second  approximation  of  z£  corresponding  to  the 
approximation  of  the  surface  pressure  by  the  first  two  terms  of  (37). 

If  this  approximation  is  considered  adequate  up  to  a  time  tQ,  then  the 
value  of  Zg  at  this  time  is  Zg^: 


*20  *  Vo  +  Vo 


(51) 


A  starting  point  z^(t)  for  the  confutations  using  Equation  (31)  is 
now  obtained  by  substituting  (51)  into  (35)  and  (36) . 
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Zt  may  be  shown  that  using  only  the  first  tern  of  (51)  i*  equiva¬ 
lent  to  the  first  approximation  given  above. 

It  should  be  noted  that  in  any  case,  in  the  successive  steps  using 
Equation  (31),  Is  confuted  first  fran  Equation  (31)  and  then  tg  is 
computed  from  Equation  (27).  The  next  step  starts  from  z^(t)  computed 
by  Equations  (35)  and  (36)  using  «g  in  place  of  a^  and  tQ. 

V.  NUMERICAL  EXAMPLE: 

The  following  data  are  chosen  as  an  illustration: 

e„  -  0.02 

c 

E  -lx  10^  psi 
p  -  k.Ok  slugs/ft^ 

From  these  data 

C  -  5»960  ft/Bec 

The  pressure  time  history  at  the  surface  is  assumed  given  as  shown 
in  Figure  6.  The  peak  pressure  pQ  at  t  -  0  is  not  shown  but  was  184,000 
psi.  The  impulse  for  the  assumed  loading  is  shown  in  Figure  7. 

The  computation  is  begun  by  applying  the  first  approximation  dis¬ 
cussed  above  for  the  first  10  milliseconds.  Equation  (31)  is  then 
applied  repeatedly  to  devolop  the  z  versus  t  curve  shown  in  Figure  8. 

Fran  Equation  (22),  £  can  be  derived  for  various  t  and  a  ,  u  can  be  com- 

S  ZB  S 

puted  from  Equations  (5)  and  (6).  Figure  8  shows  the  values  of  z  ,  c 

8  28 

and  ufl  plotted  versus  t. 

To  find  the  stress  and  velocity  at  a  field  point  z  -  500  feet,  t  - 
200  ms,  the  characteristics  through  this  point  are  drawn  to  locate  points 
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1  and  2  on  the  z  versus  t  curve  as  shown  in  Figure  8.  Reading  the  a 

3  KB 

and  u  curves  directly  below  points  1  and  2  yields  the  values 

8 


^2  -  4,000  psi 

Ug  -  59  ft/seo 

■  10,000  psi 

u ^  -  103  ft/sec 

Substituting  these  values  into  Equations  (20)  and  (21)  yields  the  field 
values  for  z  -  300  feet,  t  ■  200  ms: 

o  »  3,230  pel  u  ■  62.6  ft/sec 

The  values  of  0  and  u  at  the  surface  and  at  depths  of  500  and  700 

% 

feet  were  computed  for  various  tlmeB.  The  results  are  plotted  In  Fig¬ 
ures  6  and  9  against  time  after  arrival  of  the  Bhock  wave  for  each 
point . 

VI.  COMPARISON  OF  PLASTIC -ELASTIC  MEDIUM  TO  A  LOCKING  MEDIUM. 

The  blast  pressure  at  the  surface  a  ■  0  derived  in  the  exang>le  above 
for  a  plastic-elastic  medium  was  applied  to  a  locking  medium  also  for 
comparison  of  the  responses  of  the  two  media.  The  locking  medium  Is 
defined  by  the  stresB-strain  curve  shown  In  Figure  10. 

The  response  of  a  locking  media  to  a  given  surface  pressure  of  the 
type  shown  in  Figure  3  haB  been  studied  in  Reference  [l]  .  .  .  The  results 
may  also  be  derived  as  a  limiting  case  (c  -*  »)  of  the  problem  treated 
above.  The  location  of  the  shock  wave  front  at  any  time  is  given  by 

*.  ■  [k  / I(t>  4t  T  (52) 

o 

(1]  Salvadorl,  Mario  G.j  Skolak,  Richard,  and  Weidlinger,  Paul  "Stress 
Waves  in  Dissipative  Media"  Transactions  of  the  New  York  Academy  of 
Sciences,  Scr.  II,  Vol.  21,  No.  5,  Peps  427-434. 
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where  l(t)  is  the  Impulse  applied  at  the  surface  fjo-i  tin.-  ‘/..ro  to  time  ti 

t 

I(t)  -  J  p  (t)  dt  (53) 

o 

Figure  7  shows  the  impulse  of  the  surface  pressure  used  in  the  present 
examples.  The  stress  immediately  behind  the  shock  front  is 


SB 


P«, 


I M 

s 

.  8 


12 


(54) 


and  the  particle  velocity  for  all  particles  behind  the  vave  front  is: 


u  -  Ifcl 


Pz. 


(55) 


The  numerical  results  for  the  following  asi  umed  values  are  shown  in 
Figure  11: 


p  ■  4.04  slugs/ft^ 

e„  -  0.02 

c 

These  are  the  same  values  as  assumed  for  the  plastic-elastic  medium.  The 
vertical  stress  and  particle  velocity  immediately  behind  the  shock  wave- 
front  in  the  plastic-elastic  medium  due  to  the  same  surface  pressure  cure 
shown  in  Figure  11  also  for  ccanparison.  In  both  cases  the  stresses 
shown  are  the  peak  streBseB  experienced  at  a  given  depth.  The  same  is 
true  of  the  particle  velocities. 
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SYNOPSIS 


Asymptotic  short  tine  and  long  time  solutions  ars  derived  in  closed 
form  for  the  spherical  wave  propagation  In  an  Infinite,  Ideal  locking 
medium  due  to  various  types  of  time  varying  pressures  applied  to  the  sur¬ 
face  of  a  spherical  cavity.  Among  these  ares  (a)  a  constant  pressure; 
(b)  a  center  of  dilatation;  (o)  an  adiabatic  gas  expansion;  (d)  an  energy 
input  due  to  an  exploding  point  mass. 

INTRODUCTION 


The  propagation  of  plane  waves  in  various  types  of  locking  media  has 
been  considered  by  several  authors  w.w.w-  Spherical  waves  in 
plastic  locking  media  were  investigated  by  Konpaneets  [j  ]  and  more  recently 
by  Zvollnskil  who  derived  the  wave  equations  for  an  elasto-plastic  medium 
of  limited  compressibility  (jSj  . 

Since,  in  general,  closed  form  solutions  of  this  type  of  problem 
cannot  be  obtained,  asymptotic  closed  form  solutions  are  given  in  this 
paper  for  an  Infinite  ideal  locking  medium.  An  equation  of  motion  is  first 
derived  for  large  strains  and  displacements  when  a  variable  pressure  is 
applied  to  the  surface  of  a  spherical  cavity.  It  is  then  shown  that  its 


03  "Stress  Waves  in  Dissipative  Media"  by  M.  G.  Salvadori,  R.  Skalak, 

P.  Weidlinger,  Transactions  New  York  Academy  of  Science,  Ser.  II,  Vol.  21, 

No.  5,  pp  427-434,  1959. 

[_5 J  "  Waves  and  Shocks  in  Locking  and  Dissipative  Media"  by  M.  0.  Salvadori, 
R.  Sfcalak,  P.  Weidlinger,  Journal  Engineering  Mechanics  Division,  ASCE,  Vol. 86, 
April  i960 

[V]  "On  the  Plane  Motion  of  Sami"  by  A.  I.  Ishlinskii,  Ukr.Math.  Journal, 

Vol.  6,  No.  4,  1954. 

jjl  "Shock  Waves  in  Plastic  Compacting  Media"  by  A.  8.  Kompaneets,  Proc. 
Academy  of  Science  ,  USSR,  1956,  Vol.  106,  No.  1,  pp  49-52. 

QQ  "  "On  the  Mission  of  aa  Elastic  Wave  from  a  Spherical  Explosion  in 
the  Ground"  by  N.  V.  Zvollnskil,  P.  M.  M.  Vol.  24,  No.  1,  i960,  pp  126-133. 
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early  time  asymptotic  solution  is  tbs  rigorous  solution  of  tbs  sans  prob¬ 
lem  whan  snail  strains  and  dlsplacsasnts  ars  considered.  lbs  long  tins 
asynqptotio  solution  of  tbs  sans  aquation  la  than  sbovn  to  bs  tbs  rigorous 
solution  for  a  point  source  prsssurs  vbsn  largo  strains  and  displacements 
ars  again  consldsrsd. 

Ibs  behavior  of  a  three-dimensional  ideal  looking  nsdlua  under  com¬ 
pressive  stresses  Is  characterised  as  follows t 

(a)  At  first  tbs  material  of  Initial  density  pQ  does  not  offer  re¬ 
sistance  to  compression  but,  once  a  critical  value  of  the  strain  la  reached, 
the  density  changes  abruptly  to  a  value  and  a  constant  value  of  tbs 
dilatation: 

I-  1  -  <  1  (1) 

Is  maintained  from  then  on  (Figure  1)  . 

(b)  During  the  subsequent  motion  with  a  constant  density  p^  tbs 
material  satisfies  the  three-dimensional  yield  condition  derivable  from 
Coulomb's  lav  of  failure  [_9j  : 

ar  "  °e  "  *c  +  I  k  ^Br  +  W 

where,  due  to  symmetry,  the  radial  and  tangential  stresses  a*,,  og  are 
principal  stresses,  the  constant  c  Is  proportional  to  the  cohesion,  and 
k  la  a  function  of  tbs  coefficient  of  Internal  friction. 

tbs  above  described  behavior  Is  typical  of  certain  cohesive  granular 
soils  under  high  compressive  stresses.  Since  under  such  stresses  the  ef- 
feot  of  cohesion  becomes  negligible  as  compared  to  that  of  friction, 

GO  "On  Coulomb's  Lav  of  Failure  in  Soils",  by  R.  T.  8hield,  Journal 
Mach.  fbys.  Solids,  M1955)*  pp  10-16 ;  May  i960. 


Equation  (2)  becomes: 
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with 


e .  a  ka 

0  r 


0  <  k  <  1 


(3) 

(3a) 


2.  SPHERICAL  WAVE  EQUATIONS. 


Consider  a  spherical  cavity  of  Initial  radius  in  an  infinite 
ideal  locking  medium.  From  t  ■  0  on  the  cavity  boundary  la  subjected  to 
a  uniform  pressure  of  intensity  p;  which  is  a  function  of  the  expanding 
cavity  radius  rQ.  Let  r  be  the  Lagrangian  coordinate  of  a  particle,  and 
v  its  radial  displacement  at  a  tine  t,  when  the  medium  has  been  contacted 
to  a  density  up  be  a  dl stance  R  from  the  oavity  center,  while  beyond  R 
the  medium  la  at  rest  with  a  density  pQ  (figure  2). 

Under  these  assumptions  conservation  of  mass  requires  that: 

(r  -  w)3  -  3-— r3  -  g-j  (*0 

and,  in  particular,  that  at  the  cavity  boundary,  where  r  •  r^  +  w^-  rQ: 


e  —  -  [l  ♦  (1  -  |)(rt/h)3  ] 


•1/3 


(5) 


For  i  «  1,  and  at  an  early  time  when  the  change  in  the  cavity  radius  is 
negligible,  so  that  ra»  r^,  the  ratio  x  approaches  the  value: 

x#  -  R/r4  (*) 

At  late  times  when  R  is  large  enough  for  (1  -  t)(ri/R)3  to  he  much  less 
than  (,  the  ratio  x  approaches  the  constant  value: 

-1/3 


*L  *  I 


(7) 
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Denoting  time  derivatives  by  dots,  the  particle  velocity  at  r  (v  -  r) 
becomes,  by  Equation  (4) : 

v  -  (R/r)2  *R  (8) 

and  at  the  vave  front  (r  -  R) i 


w]  -  tR 

-r  -  R 


The  particle  acceleration,  by  Equation  (8),  becomes i 
v  -  ^(R2^  +  2RR2)r‘2  -  2l2RUR2r’5 


Tho  equation  of  motion  for  an  elementary  volume  of  soil  in  spherical 
coordinates  (Figure  3)  la  given  by: 


or  -  -  Plv 


where  v  is  defined  by  Equation  (10 ),  and  its  particular  solution  satisfy¬ 
ing  the  boundary  condition: 


rj  “  P(rc 


is  for  k  /  1/2 j 


°r  -  p^o^r) 
o 


-2+2k 


i  Pil  f?2K 

+  71 . _  i'vW  V«  R 


t  +  2RR2)  (~ 


(1  ♦  *)*: 


j  -2+2k  -1 

)  </>  •  -  (^)  + 


-2+2k 


•R  R  (-£)  -  (/) 


and  for  k  ■  1/2: 


°r  -  ptro)(r)  1_  r~  <r2* +  2R&2)(r)  + 


0  o 


2  Pl6  4*2  n  r  1  "1 

'Jj" 


(k  -  1/2)  (13a) 


•5" 


Conservation  of  noiatntun  across  tha  wave  front  (r  ■  R)  requires,  by 
Equation  (9),  that  1  r 


1(1  -  I)*8  (i^) 

Jr  -  R 

Batting  this  valua  of  oy  in  Eq'. (13)> (13*)>  and  Baking  uaa  of  Equation  (5), 
tha  aquation  govaralng  tba  notion  of  tba  vava  front  baoonaa: 


"*ir  ♦  ***  ■  ^itn) 


(15) 


In  which  tha  varlabla  ooafflciants  a  and  0  ara  glvan  for  k  /  1/2  byi 

1  /_.  .£-2kx 

*  "  Tsr^mrru  ) 

»•  .[d“. 

and  for  k  -  1/2  byt 


2  /  2-2kv  »  z  4  2 

(2k~~Tja"^'n~  (x  ‘  x  }  *  xrv'kKr'~tT  (x  *  * 


(k  /  1/2)  (16) 

-a>] 


a  *  1  -  ■ l  x  fa  x 

(k  -  1/2)  (16a) 

0  -  2  [x  ♦  x^n  x  -  |  j-f-j  (**  -  *)] 

III.  Short  Tima  Asymptotic  Solution  (Jinlts  cavity,  small  strain,  man 
displacement  solution). 

Substituting  x  ■  x|  from  Equation  (6)  Into  Equations  (16)  one  obtains 
tha  dlffarentlal  aquation  valid  for  small  displacements.  A  further  simplifi¬ 
cation  Is  obtained  by  considering  tha  strains  to  be  small,  l.e.,  neglecting 
9 

terms  In  |  and,hence,  tha  third  term  at  tha  right  hand  member  of  0  In 
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Equations  (16)  or  (l6a).  Since  Equation  {6}  is  valid  vhen  the  expansion 
of  the  cavity  is  negligible,  and  the  pressure  p  is  considered  a  function 
of  the  cavity  radius  rQ,  the  pressure  is  to  be  assumed  constant: 

P(r0)  -  pQ  (17) 


For  an  explosion  of  yield  W  and  gas  constant  J\ 

M 


p0  -  (r  -  i) 


w 


(a) 


The  solution  of  Equation  (15)  in  terms  of  the  non-dimensional  wave  front 
radius : 


x  *  R/r. 
a  '  i 


thus  becomes: 


2(2k  -  l)pr 


e  t  2 
pi£ri 


*V  -  x  1'2k)2 

s'  s  ' 


3(xb  -  *>  ' 


(6a) 


1  ~  (x^2k-  1) 


k  -  2k 


(18) 


The  corresponding  solution  for  k  ■  1/2  is; 


x  m  — 


pl«ri 


1 


x^(3/n  x  -  1)  +  1 


(18a) 


IV .  Long  Time  Asymptotic  Solution  (Point  Source,  finite  strain,  large 
displacement  solution). 

Assuming  for  x  the  value  given  by  Equation  (7),  the  coefficients 

a  and  0  of  Equation  (15)  become  constant,  and  depend  on  J  and  k  only.  This 
equation,  rigorously  valid  for  ^  -  0,  coincides  vith  the  equation  derived 
by  Kompaneets. [7] 

M  "A  Method  of  Concealing  Underground  Nuclear  Explosions"  by  A.  L.  Latter, 
•  .  .  L.I/'vicr,  E.  A.  Martinclli,  W.  G.  McMillan;  R-346,  The  RAND  Corporation, 

L'C.i  jQ }  L'0!/0* 
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■m  folXwlnc  particular  solutions  of  this  equation  are  of  interest* 
(a)  Cantor  of  BUatatlon. 

The  nagnltude  of  the  total  foeee  enerted  on  the  enrlty  ourfaeo 
la  oaaotaat  in  tlaa  so  that,  letting* 


P  -  2±>  fcvfefe? 

•  ij**!  1  1 


U9) 


tho  variable  pressure  y,  by 

t(»)  - 


k*i 


*o 

m 


(*)  on*  (9), 

_  .-a 


no  corresponding  solution  of  Equation  (19)  Is  i 


»*.!*.  - i. 


“I*  e«t2/3(»  .  a»)(i  .  i) 


»"2 


(20) 


(21) 


If  at  a  tins  t  -  whsa  R  ■  R^,  PQ  booonos  aero,  tho  wn  front  propagator 
vith  a  velocity  i 

R2  -  - 57g~ - - - -  (|-)"P^a  (R  >  Rj)  (22) 

plc  2*r/J(P  -  2a)  (l  -  0  *1  ■  * 


Blnoo  p/a  >  2  for  all  values  of  |  and  k,  R  approaches  aero  asympto¬ 
tically  and  tho  wave  front  never  stops.  This  result  stems  fron  the  neglec- 
tlon  of  the  oohealon  o  la  Rotation  (2)t  a  finite  value  of  o  gives  a  finite 
tine  at  vhlch  the  motion  stops. 

(b)  Expansion  of  An  Ideal  Gas. 

In  the  expansion  of  an  ideal  gas  the  pressure  p  follows  the  lent 


Hn  Ptr^rj7  -  p(r  )r^r  •  K  (23) 

^-0  11  00 

Shore  r  end  K  are  oonstant  characteristics  of  the  gas  and  of  the  energy 
of  the  explosion. 
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The  solution  of  Equation  (15)  vhen  the  pressure  p(rQ)  Is  defined  by 
Equation  (23)  and  p/a  >  37  ,  is  1 


2K  l-y  „-37 

p?  <T  WJ(i-i)R 


(24) 


In  this  case  also  the  introduction  of  a  finite  value  for  the  cohesion  intro¬ 
duces  a  finite  stopping  time,  as  shovn  by  Koopaneets  (jj 
(c)  Exploding  point  mass.  [12] 

Consider  a  point  mass  M  located  in  the  medium  at  r  ■  0,  which  at 
t  •  0  explodes  so  that  its  particles  acquire  a  radial  velocity  rQ.  Assum¬ 
ing  that  these  particles  remain  on  the  boundary  surface  of  an  expanding 
spherical  cavity  of  radius  rQ,  the  toted  kinetic  energy  at  a  time  t  is; 

R 

T  -  |  J helper2  dr  +  |  Mr®  (25) 

r 

0 

and,  by  Equations  (7),  (8)  and  (25),  the  rate  of  increase  of  this  energy  lsj 


dT 
dR  " 


i)r3  +MrV3 


t®R  §  +  6«p1^(f1/3  -  l)R®R® 


The  rate  of  energy  loss  across  the  shock  front  is  (Figure  4)i 


(26) 


dE 

dR 

In  the  particular  case  of  k 
medium  so  that! 


(27) 


1  no  other  energy  Iobs  can  occur  within  the 


dT  dE 
dR  dR 


O 


Hence  by  Equations  (26)  and  (27)1 


(23) 


(R3  +  rm^  +  3^  ^  “  0  (29) 

(11)  For  p/a  <  37  there  ere  no  solutions  of  the  point  source  equation. 

[12]  The  authors  are  indebted  to  Dr.  H.  H.  Bleich  for  this  solution. 
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The  solution  ef  Equation  (29)  with  R(0)  -  RQ,  and  hence  an  energy  input 
1  *9 

|  «  *0  *  iB 


where  1/2  <  4  <  1  for  0  <  (  <  !• 


VI.  Intermedlata  Range  Solutions . 

The  variable  coefficients  a  and  0  of  Equation  (15)  approach  their 

asymptotic  values  for  R — ,  but  become  practically  constant  at  finite 

values  of  R  depending  on  k  and  fc.  For  example,  the  behavior  of  a  end  0 

versus  x  «  R/r.  is  given  in  Figures  (5a)  and  (Jb)  for  k  ■  1,  1/2,  and  0 
8  X 

when  £  m  0.1. 

The  solution  of  Equation  (15)  in  the  intermediate  range  between  the 
two  asymptotic  short  tine  and  long  time  solutions  may  be  obtained  by 
numerical  means  starting  at  the  limit  of  the  range  of  validity  for  either 
solution  and  Integrating  forward  or  backward,  respectively. 

The  constant  K,  needed  to  continue  the  intermediate  time  solution  for 
an  explosion  is  obtained  by  Equation  (23).  Table  I  shows  the  range  of 
validity  of  the  two  asymptotic  solutions. 

Table  I 

(l  -  0.1) 

Short  Time  Long  Time 

5L  1.00  1.02  i.o4  1.06  1.08  1.10  1.25  1.50  10  20  - 

fi _ 

x  1.00  1.02  1.04  1.06  1.08  1.10  1.25  1.50 

8 

x  1.00  1.02  1.04  1.05  1.07  1.09  1.21  1.40 

*L . 


2.14  2.15  2.15 

2.15  2.15  2.15 


FIG.  I 
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